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TECHNICAL NOTE 3389 


AXIALLY SYMMETRIC SHAPES WITH MINIMUM WAVE DRAG 
By Max. A. Heaslet and Franklyn B. Fuller 


SUMMARY 


The external wave drag of "bodies of revolution moving at supersonic 
speeds can be expressed either in terms of the geometry of the body, or 
in terns of the body- simulating axial source distribution. For purposes 
of deriving optimum bodies under various given conditions, it is found 
that the second of the methods mentioned is the more tractable. By use 
of a quasi -cylindrical theory, that is, the boundary conditions are applied 
on the surface of a cylinder rather than on the body itself, the varia- 
tional problems of the optimum bodies having prescribed volume or caliber 
are solved. The streamwise variations of cross-sectional area and drags 
of the bodies are exhibited, and some numerical results are given. The 
solutions are found to depend upon a single parameter involving Mach 
number and the radius -length ratio of the given cylinder. Variation of 
this parameter from zero to infinity gives the spectrum of optimum bodies 
with the given condition from the slender-body result to the two- 
dimensional. The numerical results show that for increasing values of 
the parameter, the optlmimi shapes qxiickly approach the two-dimensional. 

A reciprocity relation for axial flow is derived, and it is used in 
formulating the variational problems in terms of the drag formula involv- 
ing geometry. Formulation of the minimum problems in terms of combined 
flow fields is found to lead to extremely simple relations that are sat- 
isfied by the flow fields induced by optimum bodies. The combined flow 
concepts are also useful, for example, in checking results found by other 
means . 


INTRODUCTION 


The design of minimum-drag configurations is one of iiie fundamental 
problems of aerodynamics. For many engineering purposes it is, further- 
more, possible to make useful predictions and design calciilatlons for 
steady fli^t by considering additively the drag attributable to the 
viscous nature of the air and the drag that occurs in an inviscid medium. 
Since efficient fll^t is closely associated with the use of aerodynamic 
shapes producing relatively small disturbances in the air, the analysis 
upon which the Inviscid-fluid theory is based can, in many cases of prac- 
tical interest, be further limited to first-order approximations involving 
small perturbations. For supersonic fli^t speeds such an analysis is 
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linear, the perturbation velocity potential of the flow field satisfies 
the wave equation, and the pressure drag of nonlifting configurations 
results from the accumulation of energy in the waves induced by the body 
during its motion. 

The purpose of the present paper is to show how most favorable body 
shapes, under various given conditions, can be derived by using formulae 
for drag prediction that are based upon the linearized theory. The type 
of body to be treated is a nacelle- or duct-like configuration (nonlifting 
and having axial symmetry) which Induces perturbations that are specified 
on the surface of a circular cylinder. The analysis mi gh t be termed quasi- 
cylindrical, since boundary conditions are applied on the surface of a 
cylinder rather than on the body Itself, Only the external flow is con- 
sidered. 

There are two rather different methods available for the calculation 
of drag of such bodies. The first, given by Ward in reference 1 , expresses 
the drag in terms of the geometry of the body and of a wel^ting function 
first encountered by Ll^thlll (ref, 2 ) in connection with the drag of 
fusiform bodies. The second result, published recently by Parker 
(ref, 3) > is a formula in which the -drag is expressed in terms of the 
strength of an axial source distribution that simulates the body shape. 
Generally speaking, the formula giving drag directly- in terms of geomet- 
rical characteristics would be preferable, since the usual auxiliary con- 
ditions in variation problems, such as given volume, given caliber, etc,, 
are also expressed in geometrical terms. Unfortunately, however, the 
variational problem in this case leads to an integral equation whose 
kernel is the Li^thill function mentioned previously, and the properties 
of this function are not at present well enou^ known to enable one to 
solve the integral equation by other than numerical methods. On the 
other hand, the expression for drag in terms of sources leads to. a trac- 
table integral equation, althou^ the relations between source strength 
and geometry are somewhat complex. 

Problems of the sort to be treated here have been attacked by Ferrari 
(refs, 4 and 5) by Parker (ref, 3 )« first-named author has 

approached the problem of minimum drag with assorted isoperimetrlc con- 
ditions by both the above-mentioned methods, but the main effort was made 
in connection with tiie source-strength method applied in conjunction with 
a control surface consisting of a frustum of a cone, A large number of 
cases have been worked out, mostly by nximerical methods. The other work, 
reference 3 t gives a solution to the problem of the minlmiam'-drag body 
with given caliber, making iise of boundary conditions on the Stokes’’ 
stream function, rather than the potential function. 

In this paper we shall approach the problem by the use of both methods 
outlined above. In an Introductory section, the operational approach to 
the wave equation is extended to bodies having peripheral as well as 
longitudinal variations of surface shape. The ana-lysis is then restricted 
to the case of axial symmetry and the two drag formulae are given. Then 
a reciprocity relation for axial flow is derived, and the notion of 
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combined flow fields is introduced. This device leads, throu^ applica- 
tion of the reciprocity relation and the drag formula in terms of body 
geometry, to extremely simple physical characterizations of the flow 
fields associated with optimum bodies. Next, in order to derive explicit 
expressions for some optimum bodies we consider the source-function 
approach in combination with a cylindrical control surface on which 
boundary conditions are specified. The results obtained are discussed 
with the aid of numerical examples, and, finally, the reciprocity rela- 
tions derived earlier are exhibited in terms of the explicit solutions 
found, and some uses of the reciprocity results are indicated. 

The Appendix is devoted to summarizing the results of the minimiza- 
tions for the convenience of the reader. 
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LIST OF IMPORTANT SYMBOLS 

speed of sound in free stream 
strength of source distribution 

function used in isoperimetric problems (See eqs, (60),) 

P"Po 

pressure coefficient, — 
drag 

complete elliptic integral of second kind of modulus k 
modulus of elliptic integrals 

complete elliptic integral of first kind of modulus k 
Bessel functions of order m (See ref, 8.) 
length of body 

Uo 

Mach number in the free stream, — 

S’O 

direction cosines with respect to Cartesian axes of the inward 
normal to a surface 

pressure 

pressure in the free stream 

pressure in a combined flow field, p - p 

dynamic pressure, ^ PqUo^ 
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Ar(x) 

R 

S(x) 

6S(x) 

x,y,z 

u,v,w 

Uo 

vr 

V 

Ve 

W(z) 

f 

n 

e 

n(a^,k) 

Po 

a 
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radial coordinate, •Jy^ + 

incremental radius on control cylinder due to source distritu- 
tion along axis 

radius of cylindrical control surface 
cross-sectional area of a lody 
S(x) - S(o) 

Cartesian coordiixates 

perturbation velocities in x,y,z directions, respectively 
free-stream velocity 

perturbation velocity in radial direction 
volume of body 

additional volime wrapped on cylindrical control surface 

fimctlon defined in equation (25) 

parameter of elliptic integral of third kind 

Mo® - 1 

dimensionless streamwise coordinate, y 

“X z 

angular coordinate, tan ^ 

Lagrange multipliers in Isoperimetrlc problems 

complete elliptic integral of third kind of modulus k and 
parameter a? (in notation of ref. 20) 

free-stream density 

B 

i 

perturbation velocity potential 

Suffixes 

differentiation with respect to. streamwise coordinate 
quantity evaluated in reversed flow field 
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Laplace transform 


dimensionless quantity as V* 



_S_ 


etc. 


INTRODUCTORY ANALYSIS 


The analysis to he given here is adapted to boundary conditions 
specified on a right circular cylinder so oriented that its axis is paral- 
lel to the free-stream velocity vector. Immediate application thus follows 
for quasi -cylindrical shapes that deviate sli^tly, both longitudinally 
and peripherally, from a cylindrlcaJL control surface although the expres- 
sion for drag can be extended to include the domain of slender-body theory. 

Consider a Cartesian coordinate system fixed relative to a supersonic 
free-stream velocity Uq and Mach number M© = Uo/sq > 1 where ao is 
the velocity of sound in the free stream. The x axis is alined with the 
direction of the flow and the lateral coor dinates y, z may also be expres- 
sed in polar coordinates r ,6 where r = Vy^ + 9 = tan“^z/y. A cylin- 

drical control surface of radius r = R = const, is given with the range 
0 < X < i and on this control stirface the perturbation velocity compo- 
nents, together with their gradients, are small relative to Uq and Uq/I. 
Under these conditions the field external to the cylinder of radius R 
has for its governing equation the linear relation 


P^'PxX - Tyy - 9zz = 0 (l) 

where the subscript notation denotes partial differentiation, cp(x,y,z) is 
the perturbation velocity potential yielding the perturbation velocity 
components 


>^(x,y,z) =cpx(x,y,z), v(x,y,z) = 9 y(x,y,z), w(x,y,z) = cpz(x,y,z) 

and = Mo^ - 1 . The boundary conditions on the body are to be taken 
in the form 

^r^ ^ i'eR ~ Uo^(x,6) , 0 <x < I (2) 

where G is a known function of x and 9 . 
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A General Solution of the Wave Eguation in Cylindrical Coordinates 

If equation (l) is rewritten in the form 

“ 'Prr “ = 0 (3) 

it is possible, throu^ separation of variables, to derive a general solu- 
tion representing a rectilinear distribution of source and multipole 
singularities. This general solution can be found by use of the Laplace 
transformation. By definition, the Laplace transform^ (see ref. 6) of a 
function F(x,r,S) is F(s;r,0) where 

F(s;r,e) = / e"®^F(x,r,0)dx (4) 

o 

If one employs this transformation and applies initial conditions consist- 
ent with supersonic flow theory (ref. 7) > equation (3) becomes 


3 ®s^ - - (l/r)cpj. - (l/r)^e 0 =0 (5) 

The transform of the pertubratlon velocity potential is assumed separable 
in the form 

9 (s;r,e) = ^(r,s)coB m0 

and it follows directly that ^(r,s) must satisfy the ordinary differential 
equation 

_J!L_ + J:_ + 

d(prs)® d(prs) L (prs) 




Thus, the solution can be written 


00 

cp(s;r,0) = - ^ cos ra0^i^(s)Km(PrB) + Bm(s) Inl(^^s) J 

o 

^It will be assumed throu^ the present section that the origin lies 
upstream of all disturbance points in the flow field. Subsequently, the 
origin will be shifted so as to lie at the upstream face of the control 
surface or body. 
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where Kjn and Im are modified Beesel functions in the notation of ref- 
erence 8. The asymptotic expansions for the Bessel functions show that 
yields incoming waves suitable for the analysis of flow Inside a tube 
or cylindrical control surface; Km yields outgoing waves that are suited 
to the calculation of the field external to a tube. It follows that one 
has, in the latter case, 

00 

cp(s;r,0) = - ^^Am(s)Km(prs)cos m© (6) 

o 

The inversion of equation (6) can be achieved in two ways. First, 
from reference 9^ page 277 ^ and the convolution integral, one gets 
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where the function Cin(x) is given by (from operational calculus rules) 


/ X 0X3 pXs 

dxm * ' * / <3xa / /^(xi)dx 


(m - 1)1 


i r (x - Xi)“ \a(xi)dx] 


(9) 


Equation (8) expresses the solution in the usual form, given, for example, 
in reference 10, page 527- For some purposes, numerical calculations 
for example, equation ( 7 ) has ah.vantages over equation (8). The two solu- 
tions express the perturbation velocity potential in terms of distributions 
of singularities along the central axis, the first term representing a dis- 
tribution of supersonic sources of strength. Ao(x)dx, and the subsequent 
terms representing multipoles of order m. 

It is of interest to calculate the limiting forms of equations ( 7 ) 
and (8) for large and small values of r. For large r, equation (6) 
becomes 

00 

9 (sjr, 0 ) Z - ^^^(s) J e'P^^coB m0 (lO) 

o 


where the asymptotic form 



e 


-z 


has been used. The perturbation velocity potential is then 






2 «< 72 |ir Z_r 

o 


"7x - I - pr 


( 11 ) 


The ultimate attenuation of cp with latereil distance is therefore fixed 
by the factor l/s/r. For small r, equation (6) becomes 


cp(s;r,e) « - 


2jr 


■-Ao(s) (in + /) . 


,m 


(Wi) ^ 


( 12 ) 
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vhere ' 7 = 0.5TT is Euler’s constant. Hie inversion of ecLuation (12) is 


cp(x,r,e) ~ — 


Ao(x)ln^ + ^y^ Ao(xi)ln 


X - Xi dxi + 




cos me Qm(x) 


(13) 


where Cin.(x) is defined in equation (9) . This result was used by Ward 
(ref. 11) as a basis for the development of slender-body theory. 

As presented, the above general solutions (eqs. (7) and (8)) were not 
related to specific boundary conditions. The formal development of this 
relation is strai^tforward and leads to an explicit solution for boundary 
conditions given on the cyllndrlcal-conbrol surface at r= R = const. 

Let the given conditlcJns be 


<Pr]r=R 


UoG(x,e) 


= Uo 


CO 

^gm(x)cos m© 
o 


(Ih) 


From equations (6) and (l4), one has 

00 

9r(sjr,0)]^^ = - •^^^(s) ^ Km(prs) 


r=^ 


00 

cos m© = Uo^'^(s)cos m© 
o 


Since 


^ Kbi(prs) = PsKm’Ors) 


equation (15) yields 


Am( s) — 


2jtUo lm(®) 


P sKm’(pEs) 


and the transformed velocity potential is, from equation (6), 


cp(s,r,©) - p 2. s V(3Rs) 


cos m© 


(15) 


(16) 


(IT) 
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In order to give the desired expression for cp(x,r,0) it is necessary to 
calculate the inverse Laplace transform of the functions Kni(Prs)/Kin' (pRs) . 
This task has been undertaken by Mersman (ref. 12), 


External Wave Drag of Quasi -Cylindrical Body of Revolution 
In Terms of Its Geometry or Source Distribution 


Attention is nov restricted to flow fields possessing axial symmetry 
with respect to the stream direction. Independence with respect to 9 
then reduces equations (7) and (8) to 


cp(x,r) = 


and the velocity potential is 
supersonic source potentials, 
form 


cp(s;r) = - ~ Ao(s)Ko(prs) (19) 


1 Ao(xi)dxi 


T.f 


7(x - Xi)2 - 


( 18 ) 


expressed as a rectilinear distribution of 
Operationally, equation (I 8 ) takes the 


The axes may now be considered as shifted so that the source distribution 
starts at x = - PR and induces perturbation velocities on the cylindrical 
surface r=R, O^x^Z. For r > R one then has the disturbance field 
associated with a body of revolution that deviates only slightly from the 
cylinder r = R. The wave drag of such a body can then be expressed in 
two ways: first, as a function of the body geometry; second, as a func- 

tion of the source-strength distribution. The first result has been given 
in reference 1. To the order of accuracy to which this control- surface 
theory applies, the slope of the resulting surface is 


dr ^ S'(x) 
^ 2jtR 



r=R 


( 20 ) 


where S*(x) is the streamwise derivative of local cross-sectional area 
of the body. This condition, together with equation (l9)^ yields 


Ao(s) 


Ur 


pRsKi(pRs) 


( 21 ) 
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where S*(s) means the Laplace transform of S*(x), and 


cp(sjR) _ _ S* (s) Kq(PRs) 
Uq “ 2jtPRs Ki(pRs) 


( 22 ) 


In order to calculate drag^ pressure on the body is next evaluated. 
Denoting by p and Pq local and free-stream pressure and setting 
1 2 

Qo ~ 2 linearized theory 



2u(x^R) 

Uo 


From equation (22) 


(23) 


^(sjR) 

Uo 


1 

2jtPR 


S'(s) 


Kq(PRs) 

Ki(pRs) 


S*(s) r^_ Ki(pRs) - Kp(pRs) 
2 jtPR L Ki(PRs) 


(24) 


The inverse transform of the second term involving the Bessel function 
leads to' the function W(x) introduced by Lighthill (ref, 2) , By defini- 
tion, its transform is 


_ Ki(s) - Kp(s) 
Ki(s) 


(25) 


Pressiare distribution on the body can then be calculated from the expres- 
sion 


P - Po 


<3o 


1 

JtPR 


s'(x) 

O N ^ . 

W(x) 


(26) 


The function W(x) is shown in sketch (a); 
tabular values for -2 < x < 10 are given 
in reference 1, 

The external -wave drag of -fche 

body is finally determined by direct 
integration 



Sketch (a) 
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drag 

itR^qo 



^-^S'(x)dx 


■{21) 


and from equation ( 26 ) is 


2pjf2R3 



[S»(x)] dx 


Wf f S’(i)S*(xi)w(J^-3/^*)dx dxij 


( 28 ) 


In a later section entitled “Geometric Criteria for Minimum Drag,” the 
role equation (28) plays in problems involving drag minimization will he 
discussed. For the present, it may he remarked that althou^ the magni- 
tude of the Influence function W(x) is knoTO, its analytic properties 
are not well ehou^ defined to permit easy manipulation. It will become 
more apparent later that for certain minimum-drag problems an advantage 
is provided when one deals directly with source distributions and estab- 
lishes the relationship between geometry and source strengths as a sepa- 
rate part of the analysis. 

Equation ( 18 ) expresses the potential of a source distribution of 
strength AQ(x)dx, On the cylindrical control surface r = R and within 
the range 0 < x < I an effective body shape is Induced and the drag 
of this body can be calculated as follows. The streamwlse and lateral 
perturbation-velocity components are, respectively. 


cpx(x,r) 



Ao'(xi)dxx 
•J{zc - xi)® 7 P^r® 


( 29 ) 


cpr(x,r) 



(x - xi)Aq* (xi)dxx 
7(x - xi)^ - p^r^" 


( 30 ) 


where Ao(xi) =0, for xi < - PR, The effective body, within the range 
0 < X < 1, is fixed by the boundary conditions of equation (20) and its 
external wave drag is 
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D = - 2jtPnR 


f 


cpx(x,R)(pi.(x,R)<ix 


Po ^x-3R Aq»(xj.)(x - xi)dxx 

'^ 3 R V(x - xi)2 - 3 ^r 2 


^o* (xg)(ixa 

'^PR 7(X - X2>^ - 3^^ 


(31) 


The dummy vaxlables Xi, Xg can he interchanged; if one then combines the 
two expressions of equation (3l) an<i inverts the order of integration^ the 
integration with respect to x can be performed and there results 


p po-ps^. _ 

D = — y Ao*(xi)dxi J Ao’(x2)cosh 

^-3R -3R 


2-3R 


-3R 


il - xx)(z - X 2 ) - 

3R(xi - X 2 ) 


|dX2 

( 32 ) 


as given in reference ( 3 ) . 

It is of interest to remark that althou^ equation (32) uses only a. 
knowledge of the function Ap(x) in the range -3R S x < Z - ^R^ the drag 
that is calculated presupposes a specific source distribution function in 
the range I - 3R < x if one wishes to Identify the drag with a geometric 
shape. Thus, as in sketch (b), if the body shape near r = R is assumed 



Sketch (b) 

to have some arbitrai*y variation for 0 < x < i, and to straighten out into 
a purely cylindrical surface downstream, of x = Z, a source distribution 
function is required downstream of x = Z - PR to produce the cylinder. 
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The fact that the streaM velocity is supersonic means that upstream 
influences of Aq{x) for x > Z - PR cannot he felt on the hody and 
explains why the drag of a complete geometric shape can he detemined 
from its source distribution without knowing the complete details of the 
distribution function. 


r 



As another example of the use of 
eoLuation (32) consider, as in sketch (c), 
a circular .body extending from x * - PR 
to X = Z with a cylindrical afterbody 
of radius R aft of x = Z. If the 
source distribution of this body is known 
as, say, for example, in the case of a 
cone or slender body of revolution, the 
body drag can be determined by using the 
surface r=R, 0 <x<Z as a control 
surface and calculating momentym trans- 
port through the control surface. Equa- 
tion (32) is the exact expression for 
the body drag, and, again, requires no 
knowledge of source strength beyond 
X = Z - PR, 


COMBINED FLOW FIELDS 


One method of attack that has proved to be extremely helpful in the 
analysis of problems in aerodynamic theory involves a symmetrlzation 
process in which flow fields in both forward and reverse flow are related. 
Attention, up to the present time, has been devoted principally to planair- 
type problems and in reference I3 Jones has used this approach to derive 
criteria that appear in the minimization. of wave drag of, for example, 
nonlifting wings having specified thickness ratios or volumes. In this 
section, a brief discussion is given, using the methods of reference l4, 
of the way these concepts appear in cylindrical-control-surface analysis. 


The Reciprocity Relation for Axial Plow 


Equation (l) can be written 

L(cp) s p^cp^ - cpyy - cp^z = 0 (33) 

where L(cp) is a self-adjoint linear operator. Let now ilf(x,y,z) and 
£i(x,y,z) be two solutions of equation (33) satisfying boundary conditions 
given on a circular cylinder. Reciprocal relations between and D can 
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be derived by applying Green's theo- 
rem over a prescribed geometric 
region. Consider, as shown in 
sketch (d), the cylindrical control 
surface extending from x = 0 to 
X = 2 and draw the enveloping Mach 
cones at the front and rear of the 
surface. Denote the cylindrical sur- 
face by the fi*ont Mach cone 

X - pr = - 3E by Zg, and the rear 
cone X + pr » 2 + PR by S3. These 
surfaces enclose a toroidal region, 
bounded internally by Z^ and exter- 
nally by Zg and Z3. It follows from 
Green's theorem that the integral 
relation 

Sketch (d) 

// i ^ I * 




applies where the surface Integration extends over Z^, Zg^ Zg ytnfl m, 
na, ng are direction cosines, with respect to the x, y, z axes, of the 
surface normal directed inward into the region. 

It is customary to re-express relations like equation ( 34 ) in terms 
of a newly defined directional derivative along a line termed the conormal. 
In this manner, the equation becomes 




(35) 


where 


l 6 
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and the direction cosines Vi, Vs, V 3 of the conormal are derived from 
-niP^ = Avi, ns = Av2 , ns = Avs 

By calculation of the respective normals ni, ns^ using the relation 

Vj_2 + + Vg2 _ it is readily found from the equations defining the 

conormal that on the, surface 2^^ the conormal is normal to the surface 
and A = 1; on a Mach cone, the conormal lies along the cone and A = 3 . 

Let now t he set equal to cp(x,r,e) , the perturbation velocity 
potential, associated with boundary conditions in a forward-flowing stream, 
and let ft be u(x,r,e), the x-wise component of perturbation velocity 
associated with boundary conditions in a stream flowing in the reverse 
direction. Under these conditions, equation (35) becomes 


y's |e ax + p yja ^ //“ I? '^3 

o o 


“ /''f i ^ a//i> i ^ i “3 

o o 


On the Mach cone Eg, the perturbation potential may arbitrarily be set 
equal to zero and its conormal derivative along the cone will also be zero; 
as a consequence, the second terms on both sides of the equation vanish. 
Since the flow fields are Irrotational, Su/Sr = 5vj./8x where Vj. is 
radial velocity. After making this substitution and integrating the first 
term in the ri^t member by parts , one gets 


p 2 Jt pZ 


r \ r^~ 

r de r 

uvi-dx = R / 

cp( Z,R,e)vr( Z,R,e) - / VrU dx 

j j 

0 0 


J 

0 


o X=Z/2 

The last integral becomes 



_d^ 

dv 



dv 
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and for tlie given boundary conations it is possible to show that along a 
conormal of Eg the relation Vj* = Pu holds and sfr u is Independent of 
V. The Integral then can be rewritten, as 


o x=l/z 


PR 


J ' 


u(i,R,e)cp(z,R,e)d0 


and one has, finally, the desired reciprocal theorem 


-R 


r-f 


I 

u(x,R,e)vr(x,R, 




x,R,©)vr(x,R,e)dx 


( 37 ) 


It is not the purpose here to exploit the various applications of 
equation (37)j rather, the role played by the reciprocal relation in drag 
caJ-culations will be considered. In the forward and reverse flow fields, 
the pressTire-velocity relations of linearized theory are 


P - Po = “Po^o^^ P - Po = PoUqT^ (38) 

If, furthermore, thickness distributions of the form 


r = f(x,e), r = f(x,0) 


are placed on the cylinder r = R, the boundary conditions are 


Uo 5rJ 


r=R 


ax " 


Equation (38) can then be written 


1 a? 
^ arJ 


'r=R 


ax 


p23f pi ^ p25C pi 

J R de J (p - p^) |- d^ . J R dej (p 


( 39 ) 
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An immediate consequence of this last result is that for a unique 
thickness distribution, that is, for f(x,0) = f(x,6), the drag of a 
body in forward and reverse flow is the same. This follows from the fact 
that for quasi -cylindrical bodies the relations for drag are, respectively. 


D 



PoJr=R 


dx 


o o 


nZTt pi 

E = -eJ isj 

o o 

For fixed geometry, therefore, drag is equal to. half the sum of these two 
' expressions 


Defining pressure P(x,r,6) in the combined flow fields by the following 


P(x,r,e) = p - p = -PoUq(u + u) 


.M. 


one has 




P(x,R,0) dx 

OX 


m 


If the body has axial symmetry, equation (^l) reduces to the form 
given in equation (28), To show this, one notes first that P and f are 
independent of 0 and that eq-uation (4l) becomes 
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The proof follows from the relations 


S»(x) = S»(x) 




^(x,R) 

Uo 






> (h2) 


P(x,R) ^ 
9o 






Geometric Criteria for Minimum Drag 


Consider now the problem of minimizing the wave drag of a quasi- 
cylindrical body subject to the condition that the volume of the body 
is constant. The body surface may be defined by 

r = f(x,e) = R + g(x,e) (I1.3) 

The function g(x, 0 ) determines the magnitude of the surface displacement 
from the cylinder r = R; these displacements, as well as their gradients, 
are assumed small and we also assume 


g(x,e) = 0 for X < 0 and Z ^ x 


If equation (If-l) is integrated by parts with respect to x, the wave drag 
of the body becomes 


D 



I 

P*(x,R,0)g(x,e)dx 


m 


where the prime indicates x-wise differentiation of P. The imposed 
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geometric constraint on the variational problem is 


I r-i 


f^(x,0)dx 


i r-f 


[R2 + 2Rg(x,e)]dx 


= nR^2 + R 




I 

g(x,e)dx 


= V 


const. ( 45 ) 


where V is the total volume. The problem thus becomes one of minimizing 
the expression 


D - liV = 




I 

P‘g(x,0)dx 


+ 





where u is- the Lagrangian multiplier. Carrying out the variation, one 
has 


5(D - pV) 



+ 2p6g]dx = 0 


but from equation (37) or (39) it can be shown that the first two terms 
in the integrand yield equed. integrals and the minimizing condition 
becbmes 



I 

[P»(x,R,0) + ti]6g 


dx = 0 


Since this latter equation must be satisfied by all possible varia- 
tions of the displacement function g(x,0), it follows that the desired 
condition is 

p*(x,R,e) + p = 0 (47) 

stated in words, the condition for minimum wave drag of a quasi -cylindrical 
body of given volume is that the longitudinal gradient of pressure on the ' 
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body in tbe combined forward and reverse flow field is a constant. Fur- 
thermore, from equation , minimum drag is then given by 

%in = ^ (V - ^ Ve (48) 

■vdiere [i can now be identified with the negative pressure gradient in 
the combined field and (V - = Ve is the volume exposed to the fluid 

aro^lnd the cylindrical control surface r = R. 

The actual cross-section-area variation of a minimum-drag constant- 
volume body and its pressure distribution are shown in sketch (e) for the 



case in which axial symmetry is imposed. Pressure coefficient in the 
combined flow field of an axially symmetric body has been given in equa- 
tion (4-2) . The geometric criterion Just established then leads to the 
Integral equation 


2S*(x) 



-px + b 


(49) 


and the solution of this equation will determine the body geometry. In 
the following section an analogous integral equation will be derived but 
with the source -strength distribution chosen as the fundamental dependent 
variable . 

The combined-flow-field technique can also be used to study the 
problem of minimizing wave drag for specified body caliber or, more 
generally, when the body has a fixed cross section at a specified 
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longitudinal position. The resulting condition for minimum drag is that 
the pressure distribution on the body in the combined flow field is a 
constant forward and aft of the specified position. These conditions 
are all analogous to those obtained for planar problems by R. T. Jones 
(ref. 13). 


DRAG MIinMrZATIOR 


In this division, optimum bodies having certain prescribed geometric 
properties will be determined by standard variational methods . The ansS.- 
ysis will, as mentioned previously, deal with” the strength of an axial 
source distribution as the minimizing function, rather than the geometric 
quantity, cross-sectional area. Thus, we shall be concerned with formula 
(32), giving drag in temns of the source distribution. 


Quasl-Cyllndrical Body of Revolution of Given Volume 


Isoperimetric conditions . - The configuration to be considered, 
together with associated nomenclatvire , is shown in sketch (f) . The geo- 
metric properties of the body 
can be expressed in terms of 
the sovirce distribution func- 
tion Ao(x) by using eqoia- 
tion (20) , namely. 



-► X 


s‘(x) = 

Then, from equation (30) 


(50) 




r^-PR Aq«(x3.Ux - Xi) 

-pR n/(x - Xi)^ - P^R*" 


dxi 


(51) 


If equation (5I) is integrated x-wise it is seen that 

,z-pR^ 


/ V / ^ 1 r^"pR(z - Xx)Aq*(x 

i(x) - S(o) = i / dz / \ 


(x i)dxi 


(52a) 
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By changing the order of integration and performing the integration with 
respect to z one finds 


S(x) - S(o) = i 
Uq 



(x - Xi)Ao(xx) 
/(x - Xi)2 - 


dxx 


or, integrating hy parts. 


(52b) 


S(x) - S(o) = ^ J Ao’(xx)^/(x - xx)^ - dxx 

-m 


(52c) 


The magnitude of the additional volume wrapped around the cylinder is 




[S(x) - S(o)]dx 


(53a) 


and, from eqioatlon ( 52 b) , 


^ J' ■Ao(^i)*'/(^ “ dxx 

° -PR 


(53b) 


The variational problem .- The quantity to be minimized will be tahen 
as D - [J.xVe. From equation (32), the drag can be written, after an 
Integration by parts. 


D 


_pjo r 

~ hjc J 


Ao(xx)dxx 


-PR 


- xi)"" - P^"" ^ 


ol-pR Ao»(x2)7[z~- xa)^ - P^^ 

J 


Xx - X 2 


pR 


(3k) 


In addition to prescribing the volume added to the fundamental cylinder, 
we shall also require that the body return at the end to the same cross- 
sectional area as at the .front. Thus, according to equation (52b), 


Z-PR (i _ xx)Ao(xx)ixx 
^/(^ - xx)^ - p^R^ 


0 


(55) 
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is a condition to be met by the minimizing function Ao(xi), and by its 
variations . 

The quantity D - iJiiVe can be formed from equations ( 53 ) and ( 5 ^)# 
and if the variation is perfomied, one finds the condition 


^-^(xx)dxi 
7(1 - xi)2 - P^r2 


Ao*(3{e)7(l - Xs)® - 
/ ~ I ^ < 3^2 


ji[U - xi)® - p^R®] 

Po^o 

. 


= 0 


If this last equation is compared with equation ( 55 )# I't is seen that for 
admissible variations, the quantity within the brackets must be set equal 
to A(z - xx) , where A is an arbitrary constant. Thus, the equation for 
determination of the optimizing source distribution under the conditions 
of given volume and closure is 


Aq'(xi) 7 (Z - xx)^ - 

X - Xx 

-pR 


dxx = A(j 


x) 


PoUo 


1(1 - x )2 - P% 2 ] 

( 56 ) 


Equation (56) is recognized as the familiar airfoil equation with 

[ Aq ’ ( xx) 7 (I - xx)^ - P^R^] as the unknown. Thus we write the solution 
immediately as (see, e.g., ref. 15) 


Ao *^ ( xx)*^ '( I -x) ^ - P^R^ = 


It 7 (i-pr-x)(x+pr) 


vl-pR 


It J' Ao’ (xx)7( l"Xx)^-P^R^'ixx- 


-pR 



A(l-xx)+^[(l-xx)^-P^R^l 

1 -PR-Xx) (xx+pR) dxx 


The first integral on the rl^t vanishes according to the closiire condi- 
tion (55) and, if the remaining integrations are performed, we find 
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Ao'(x)V(Z-x)"-P"r" = ^ + 


l-QR-x 

8 


( I -2x)+^^ ( 31 ^ -4pR2 -83^R^ -121 x+8x^ ) 
Po^o 


} 


It will be noted that unless 


2]tZ 

PoUo 




+ 


2X = 0 


this solution for Aq’ ( x) does not obey the closure requirement. There- 
fore we impose this last condition and finally obtain 


Ao*(x) 


^oUo 


(l^ - UpRl - 8p^R^) - 8Zx + 83{^ 
V(2 + gR - x)(x + PR) 


( 57 ) 


The strength of the minimizing source distribution Aq(x) is now obtained 
by integrating equation (57)j 


Aq(x) = {I - 2 x) 7 (rTir^nj( 7 TpS^ - 232 r 2 cob “1 ^ 

2poUo L I + 2 $rJ 

(58a) 

Properties of the optimal source distribution .- It is convenient to 
express the various quantities such as source strength, area distribution, 
etc., as dimensionless functions of the dimensionless variable tj = x/i 
and parameter a = pR/z. Thus, indicating a dimensionless function by a 
star, we have from equation (58a) 


Ao*(il) = (1 - + ff)(l - Ti + ff) 

L Uq L 


2o^cos*^ 


1 - 2n 


1 + 2ff J 

(58b) 
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It will be noted that if the radius of the control surface is taken very- 
small, so that CT 0 , formula ( 58 b) becomes 


Ao*(n)]o-s^o = - 2n)7n(i - n) 

which is the well-known slender-body theory result for -fche source distri- 
bution corresponding to an optimum body of given volume (refs. I 6 arid 17 .) 

In order to determine -the value of -the Lagrange multiplier (J-x, in 
terms of the prescribed volume Ve# it is convenient to find first the 
expression for the local cross-section area of the optimum body. Thus, 
using equation (52c) (wi-th S*(ti) = ~ S(Zti)), 

z 


S*(n) - S*(o) = ^7(n + 2cr)(l - n + 2ct)[ti(1-ti)E - cr(l - 4cr)(K - E) ] 

' (59) 

where K and E are elliptic integrals of the first and second kinds, 
respectively, of modulus 


k2 = 


t( 1 - Tj) 


(n + 2cr)(l - T1 + 2cr) 

Using equation (59) 1°- equation (53a) ^ we find 


Ve* = ^ f .7(11 + 2d) (l - T] + 2cr)[n(l - n)E - o(l - 4d) (K - E)]^^ 


= ^ B(cr) ( 60 a) 

which expresses "the constant p-x^ in terms of the prescribed volume 

Ve* and of a function B of -the quantity a = PR/I. A graph of 

this function B(c) versus o is shown in figure 1. Sho-wn also in 
figure 1 is a dashed line that corresponds to the asymptotic form for 
B(a), which is 

B(») « H 

The closeness of the asymptotic values to the exact values even for 
relatively small values of a is noteworthy. 
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The formulae (58b) and ( 59 ) for the source strength and cross-section 
area, respectively, can now he recast in terms of prescribed quantities 


Ao*(n) = 2 B (g) “ 2 ti)7(ti + 0)(1 - q + a) - 2o^cos’ 


-i l - 2 t ) 
1 + 2cr 


(61) 


^ +2tr)(l- n +2a)[Ti(l- ti)E -or(l -1 ut)(K- E) ] 

(62a) 

where Vp is the volume of the original cylinder section, 

Vq = jtR^Z = ZS(o) = 2 ®Vq* 

Consider the expression (61) for the soiarce function Ao*(ti). In the 
parameter a = pR/i, we may think of p as fixed and I as unity, so 
that variations in a amount to variations in the size of the control 
surface of radius R. Thus, in sketch (g), the case a = 0 corresponds 



to the source distribution for the well-known Sears-Haack body (refs. I6 
and 17) . It will be noted in the cases where a > 0 that the source 
functions become less steep and attain lesser maximum values because the 




28 


KACA TN 3389 


volume remains the same while the control- surface cylinder is increasing, 
thus giving a smaller maximum radius of the added portion. 

Next let us examine the expression (62a) for cross-sectional area. 
First, we notice that it can he written 

y * 

S*(ti)-S*(o) = J(r\ + 2a)(l- Ti+2(T)[Ti(l-T))E-cy(l- 1 kt)(K-E)] (62b) 


in which form it reduces formally for a 0 to 


S*(t)) = ^ Ve*[t)(l _- (62c) 

which is identical with the expression for cross-section area of a slender 
optimum body of prescribed volume (Sears-Haack body) . Of course, Vg is, 
in this case, the total volume of the body. On the other hand, if we 
allow the radius of the control surface to increase indefinitely , equa- 
tion (62b) gives (using the asymptotic form for B(ct), eq. (60b)5 


S*(ti) - S*(o) = 6Ve*n(l - r\) 

In the case when R is very large, we take 

S(x) - S(o) = 2JtRAr(x) 

so we have, returning to the original variables, 

Ve x(Z - x) 


Ar(x) = 6 


2rau 


(63) 


( 64 ) 


where Ve/2jtRZ is a finite quantity, and. In fact, is the average height 
of the protuberance above the control cylinder. This result is clear 
from physical reasoning, for one would expect that as the control cylinder 
increased in radius, the two-dimensional result for the optimum problem 
would become more nearly valid, and, indeed, equation (64) is the formula 
for a two-dimensional biconvex section, where Ar is distance from the 

~ • A Ve' 

mean line, I is chord length, and maximum thickness is ' 




It will be. noted that the area distribution as given by equation 
(62b) has fore-and-aft symmetry, since the functional dependence upon t] 
Involves only the combination t)(1 - Tj) . The maximum cross-section of 
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the optlmvcm hody then occurs at the midpoint ti = 1/2 and is given by 
(from eq. (62b)) 


- S*(o) = 2 V, 


1 + 4cr 
4 B(a) 


"E 

k " 


g (1 - 4 a) (K - E) 


= 2 Ve*T(a) (65) 


where the modulus of the elliptic integrals is now k = l/(l + 4 a). Sketch 
(h) shows the function T(a) versus a. 



The drag of the optimum bodies can now be evaluated. From equations 
( 54 ) and (56) 


D = 


Po Ao(|)d| 




/ 

-PR 


V(J - 


m - 






The Integral involving A vanishes because of the closure condition. The 
remaining integration gives 

.I-PR , 

D = ^ r Ao(|)7(2 - 1 )^ - ( 66 ) 

°’^R 

by equation ( 531 i) . Finally, using the evaluation of of equation 

(60a), we have 
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i = 3 

Qo l*B{a) 


(67) 


Numerical results pertaining to the problem Just solved will be given 
in a later section, and a summary of the important formulae is given in 
the Appendix. 


Quasi-Cylindrical Body of Revolution With Given Caliber 


The variational problem .- For this problem, we prescribe the area 
at the base of the body, so the given condition is, from equation (52b) 


AS = S(l) 


S(o) 



(Z - xx)Ao(xi)dxi 

sl{l - P^R^ 


(68) 


The variation can be taken an before (now without invoking the closure 
condition) on the quantity D + MS, and it leads to the integral equation 


Ao*(xi)V(Z - xi)^ - 0%^ 

J X - XI 

-m 


Po^o 


{I 


x) 


(69) 


The solution to equation (69) consistent with the given conditions is 


4 Uq(^) Z - 2x 

^ Z(Z + UpR) x/(z + PR - x) (x + PR) 


(TOa) 


Integrating this expression, we find for tb-e strength of the optimizing 
source distribution 


Ao(x) = ^ _ U o (^ S) ^(1 + PR - x)(x + PR) (TOb) 

^ Z(Z + iipR) 


The source distribution of equation (TOT^) represents the first approxima' 
tion to the result, of reference 3 for nearly equal front and rear radii. 


Properties of the solution .- As in the section on the body with pre- 
scribed volume, we now consider x made’ dimensionless by division by Z, 
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and again set 0 = pR/l. The various quantities of interest in connection 
with the caliber problem then become 

= I + 0)(1 - n + 0) ( 70 c) 


S*(ti)-S*(o) 

As* 


2 

3 t(l+ 40 ) 


(1+20) (l+lw)n(g^,k)-(l- 2 ^) (n+2g) (l-Ti+20)E-(l+20) (1 -ti+20)K 

7(ti+20)(1-ti+20) 

C71) 


where II(a? ,k) is a complete elliptic integral of third kind of modulus 
k^ = r parameter a® = — — . Again K and E are 

(ti+ 20) (l-Ti+20) 1-T]+20 

complete elliptic Integrals of the first and second kinds, respectively, 
of the same modulus k. 

If we allow 0 to approach zero, equation (?l) becomes, in the limit, 

= I |^sin“Yn - (1 - 2 ti)^ti( 1 - q) j (72) 

which is the shape function for the well-known Karm^ ogive (ref. I8) . 

At the other limit, when 0 — > 00, equation (7I) gives (in the original 
variables) 


S(x) - 3(0) _ X 
AS I 


or, using the approximation of equation (63) , 


Ar(x) _ X 
Ar( 1) I 


( 73 ) 


which is again the expected two-dimensional result for specified caliber. 

The drag can be found by substituting equation (69) in equation ( 5 ^)^ 
and then using equation ( 52 b) . There results 
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D 1^- (^S)^ 

“ it(l + Ita) 


( 7 ^) 


A summary of formulae pertaining to this hody will be found in the 
Appendix. 


Examples of Optimum Bodies 


The optimum body of given volume .- In order to examine in detail the 
dependence of the body geometry on the parameter o, we may return to 
equation (62a). The quantity [S*(tj) - S*(o)]' is actually the local cross- 
sectional area added to the basic cylinder~by the action of the soinrce 
distribution. In figure 2 are shown some cases of optimum bodies, having 
equal additional volimie Vg*, for several values of the parameter cr. 

Only half of each distribution is shown, since they are symmetric about 
the point tj = l/2. The one labeled 0 = 0 is the Sears-Haack optimum 
body, and it will be noted that as 0 increases, the curves depart rather 
quickly from this limiting case and approach the other limiting value of 
the biconvex distribution for 0 -^00. in fact, a biconvex arc drawn 
throu^ the end points of the 0 = l/2 case is indistinguishable from the 
exact result in the scale used. In the inset of figinre 2 is shown the 
variation of the drag of the optimum bodies as a function of 0. This 
drag is also based on equal volume, and shows a fairly rapid decrease 
with increasing values of 0, due to the decrease in the thickness of 
the exposed portion of the body. The dashed curve on the drag plot is 


the calculated drag 


D _ 

<lo(W^^)^ itcrj 


under the assumption that each 


meridian section of the body acts as an independent two-dimensional opti- 
mum airfoil. This admittedly crude approximation is of course very poor 
at low values of 0, but its accuracy becomes surprisingly good for 0 
greater than about 0,4, and the approximation becomes exact in the limit 
0 — > 00. 


The variation of local cross section with o can be examined also 
on the basis of equal exposed area. Thus, using equation (65) in combina- 
tion with equation (62b), we have 


AS*(ti) 

'^*max 


1 

2T(0)b(0) 


V(tj + 20)(1 - ti + 20)[ti(1 - ti)E - 0(1 


- lw)(K - E)] 
(75) 


Figure 3 shows plots of equation (75) ^ and it is again noted that the 
departure from the slender-body approximation (0=0) is rapid. The 
limiting variation of area for 0 — > w is also shown in figure 3j an^ ii 
is seen again how closely the optimum body-shape functions approach this 
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limiting result even for moderate values of o. Also shown is the drag 
corresponding to these cases. 


D 3 

c[Ql^(ASmax)^ Iffl(cr) [T(cr) 


(76) 


which shows a similar drop from the 0=0 value as 0 is Increased. 
Agaln^the effective fineness ratio of the bodies is increasing with a, 
and, if frontal area exposed to the stream is held fixed, the maximum 
thickness of the excrescence vanishes as I/0 for large 0. The departure 
of the geometric variation from the slender-hody case is most pronounced 
near the nose, T| = 0, where the slope is given by 

^ (1 + 1^) /zm (77) 

dlj 8 B(0) V 1 + 20 

which vanishes only as nTc for 0 — >0. 

The optlmtmi body of given caliber . -■ In this case, the maximum cross- 
section occurs at tl = 1 so there is no longitudinal symmetry. Figure 4 
shows, for several values of the parameter 0, the optimum, equal -caliber, 
incremental cross-section area given by equation (71) • The inset shows 
the drag as a function of a; from equation (7^) 


<lo 



h 

jt(l + 40) 


Again in this case, the closeness of the optimum distributions as cr 
Increases to the two-dimensional value (0 — >-oo) is noticeable. This 
point has also been made by Ferrari in reference 5 where problems similar 
to ours are treated by a different approach. If the expression for cross- 
section area (eq. (71)) is expanded in powers of I/0, it is found that 

AS*(n .0) ^ T](1 - T])(1 - 2 n) 

AS*( 1 , 0 ) ^ ■ 32 

which shows the smallness of the correction to the two-dimensional result 
for moderate values of 0. 




34 


MCA TN 3389 


Reciprocity Relations 


The optimum body of given volume .- The longitudinal and radial per- 
turbation velocities can be determined by substituting the derivative of 
the source-distribution function (eq. (57)) into the formulae (29). We 
find, at any point (x,r) (r > R) in the field 


9x(x,r) 

Uo 


Site 


4[x+p(r+R) ] [Z-x+p(r+R) ]E + 


si [ x+p(r+R) ] [ I -x+3(r+R)l 
[ I ( Z+4pR) -4( Z+23R) ( I -x+pr+pR) ]K+4( l+2pE) ( l -2x)II(a® ,k) J- ( 79a) 


_ 1^1 


Uq SjcOor ^[x+p(r+R) ] [ l-x+3(r+R) ] 


{(. 


2x) [x+p(r+R) ] [Z-x+p(r+R) ]E - 


PrZ[Z-2x+2p(r-R)]K-4p^R(r-R)[l-x+3(r+R)]K + 

4p2(r2-Ra) (i+23R)n(a2,k)|- (79b) 

where now 


^ x(Z - x) 

[x -jr p(r + R)][Z - X + p(r + R) J 


X - 3(r - R) 

I - X + p(r + R) 


For the present axis system, the act _of reversing the flow amounts 
to substituting I - x for x, and, for ttie case of the symmetric body, 
the longitudinal perturbation velocity in the reversed flow is 


u(x,r) = -u(l - x,r) = - " ^^^) 


How, from equation (4o), pressure in the combined field is given by 

P = -PoUo(^ + u) = -PoUo[q>x(x^r) - ~ x,r) ] (80) 
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Substituting equation (79a) into this relation, we. find 

P = (2 - 2x) ( 81 ) 

by xising the addition formula for the elliptic integral of third kind. 
Differentiating equation ( 8 I) we find 

P» + 111 = 0 ( 82 ) 


which agrees with the criterion for minimum drag with given volume estab- 
lished in equation (^7) . The Lagrange multiplier 11 1 is therefore 
identified as the pressure gradient in the combined flow field. It will 
be noted that equations ( 81 ) and ( 82 ) hold everywhere within the envelop- 
ing forward and rearward Mach cones of the quasi- cylindrical body (see 
sketch (d) ) . 

Now considering the radial component of perturbation velocity cpy, 
we find 


q)j.(x,r) qpj.(Z - x,r) - R^) 

+ Hi (83 

Uo Uq llq^r 

so that the relation 

cp^(x,r) = -cpr(^ - x,r) 

is satisfied on the quasi-cylinder itself, that is, when we set r = R. 

The opt^TnuTn body of given caliber .- For this case we find the fol- 
lowing equations for the perturbation velocites: 


Uo 


Ul + 2SR) 


2Kqa*Ax + p(r + R) ][2 - x + p(r + R) ] 


[K - 2n(a2,k) ] (84a) 


~ 


2Jtqo^tx+p(r+R) ] [ Z-x+P(im-R) ] 



][Z-x+| 3 (r+R)]E - 


pr(Z+2pR)Kl 


J 


(84b) 


where the elliptic integrals . have the same modulus and parameter as in the 
previous section. 
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In this case_, the pressure in the comhined flow field is 


P = “ x,r)l (85) 

which gives 

P = A (86) 

so that in this instance, pressure itself is constant in the combined . 
flow field. ’ - . 

From equation (84b), we see that 

9j.(x,r) = - x,r) (87) 

since the modulus of the elliptic integrals is invariant to the change 
X — > Z - X. 

Uses of the reciprocity relations .- The reciprocity relations serve 
the dual function of checking the derived perturbation potential against 
minimization criteria based on other considerations (see eq. (47)) and of 
relating the Lagrangian multipliers to the pressure or pressure gradient 
in the combined flow field. Equations (8I) and (86) also reveal that the 
expressions for pressure in the combined flow field hold. Independently 
of r, throughout the entire region within the enveloping cones of the 
bodies. These results are generalizations of a similar effect noted in 
reference 19^ where the combined pressure*-field associated with a Sears- 
Haack body was shown to have a constant gradient within the enveloping 
cones. In the latter reference, this property of the minimum-drag body 
was used to expedite the calculation of inteference drag with a satellite 
body lying within the enveloping cones. Similar methods could obviously 
be applied to the present configtirations. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Nov. 22, 1954 
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APPENDIX A 


SUMMARY OF FORMULAE FOR IHE OPTIMUM BODIES 



The formulae derived In 
the text for the body shape func- 
tion, pressure coefficient, and 
drag of optimum bodies having 
given volume or given caliber 
are repeated here for conven- 
ience. The type of configura- 
tion treated, and the nomencla- 
ture , are shown in sketch (Aa) . 


The Optimum Body of Given Volume 


The variation of AS for the optimum body with given volume is 


AS(x) 


where 



*/(x+2pR)(Z-x+2PR) -|^x(l-x)E(k) - 3R(2-1^PR) [K(k) -E(k) ]j- 


(Al) 


AS(x) = 
Ve 

<f)- 

K(k) = 
E(k) = 
k2 


«[ (R + Af)^ - R^] 
volume of exposed portion 

B(ct) function defined in eqviatlon (60a) and shown in figure 1 
complete elliptic integral of first kind of modulus k 

complete elliptic integral of second kind of modulus k 
x(Z - x) 

(x + 2PR)(l - X + 2PR) 


Examples of optimum bodies for a few values of the parameter BR/z 
shown in figures 2 and 3 . ' 


are 
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Cp = 


The pressure coefficient on the body Is 
P - Pq 
<10 



4 (x+2PR) ( I -x+2PR)E-[ 1 ( Z+4pR) -k{ Z-fgpR) ( l-x+2pR) ]K+U( Z+2fiR) ( Z-2x)Il(a^,k) 


«7(x + ^R) { I - X + 23R) 


(A2) 


where n(oc?>k) is a complete elliptic Integral of third kind of modulus 
k and parameter a® (in the notation of ref. 20) . The parameter a® 

Is given by 


_2 _ X 

Ct = - ■ - ■■■ 

I -X+2PR 



x/l 

Sketch (Ab) 
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Ttie wave drag of this optimum body is given by 


-2- = 3 —5^ 

’ i^(f 


(A3) 


The variation of drag with pR/Z is shown in figures 2 and 3* 


The Optimum Body of Given Caliber 


The veiriation of ZiS for the optimum body of given caliber is 


AS(x) 


2 AS(Z) ^ 
* 1(Z+4PR) 


I ( Z+2PR) ( Z+4pR)n(g^ ,k) -( Z-2x) (x+2pR) ( l-x+2pR)E-l ( Z+2pR) ( Z-x+2pR)K 

•/(x + 2^) ( Z - X + 2PR) 


(.Ah) 


where the symbols have been defined above. Examples of optimum bodies for 
a few values of the parameter pR/z are shown in figure h. 

The pressure coefficient on the body is given by 



AS(Z) Z + 2PR 2n(g^,k) - K 


Z + l^PR 's/(x + 2PR) (Z -X + 2PR) 


(A5) 


Z 
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